Correlations between cosmic strings and extra relativistic species 
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The recent observation that the Cosmic Microwave Background (CMB) prefers a neutrino excess 
has triggered a number of works studying this possibility. The effect obtained by the non-interacting 
massless neutrino excess could be mimicked by some extra radiation component in the early uni- 
verse, such as a cosmological gravitational wave background. Prompted by the fact that a possible 
candidate to source those gravitational waves would be cosmic strings, we perform a parameter 
fitting study with models which considers both cosmic strings and the effective number of neutrinos 
as free parameters, using CMB and non-CMB data. The implications are twofold: on the one hand 
cosmic strings may be the extra source of the gravity wave background needed to fit the data; and 
on the other, due to correlations between parameters, a lower extra radiation component may be 
needed. We find that there is in fact a correlation between cosmic strings and the number of extra 
relativistic species, and that strings account for at least a part of the extra radiation necessary, but 
it depends strongly on the cosmological data used. In fact, CMB data prefer strings at a 2a level, 
paying the price of a higher extra radiation component. CMB data also give a moderate preference 
for a model with n s = 1. The inclusion of non-CMB data lowers both the preference for strings and 
for the extra relativistic species. 
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I. INTRODUCTION 

The analysis of Cosmic Microwave Background (CMB) 
has been invaluable for our understanding of our uni- 
verse. There is now a standard cosmological model that 
the community agrees upon, which fits the data very well. 
It is in the link between the cosmological standard model 
and the high energy standard model where most of the 
unanswered questions lie. It is in the small details and 
small windows that the data allow us to try different ap- 
proaches where a promising road to understanding the 
constituents of the cosmological standard model lays. 

As an example of what we mean we can consider cosmic 
strings P, Q. These are interesting entities that many 
of the possible high energy physics models for inflation 
predict to have been formed in the early universe, and 
their detection could give invaluable information about 
the physics of the early universe. The CMB anisotropy 
coming from strings has been constraint to be less than 
10% |3-lig|. and predictions on CMB signals at small an- 
gular scales [ll|, [IH and on B mode polarization fl3l — 
161 ] show that new CMB measurements will further con- 
straint models with strings. 

Another example could be the observation that several 
recent analysis on cosmological data seem to indicate: it 
has been showed [13, Ell that the effective number of 
neutrinos N e f{ is more than what one expects if there 
are only three different neutrino species. This could be 
interpreted as more than three neutrino species being 
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present in the early universe, which seems to contradict 
what we know from high energy physics. But it can also 
be understood as the effect of having some extra radiation 
component in the early universe. 

Along these lines, some groups have been investigating 
the option of gravity waves being the extra component 
[l9| . Those gravity waves could account for the extra ra- 
diation component needed to explain the observed N e ff. 
The Cosmic Gravity Wave Background (CGWB) behaves 
as non-interacting massless particles [20j . so their effect 
on the CMB is the same as massless neutrinos when 
their energy-density perturbations are adiabatic, as one 
would expect from inflation. In other scenarios, where 
the perturbations on the CGWB energy-density are non- 
adiabatic the effects on the CMB from CGWB will differ 
from those produced by massless neutrinos [2l[ ■ 

One important question, though, is to give the source 
of those gravity waves. If they were there at the be- 
ginning, what was creating them? One of the possible 
answers put forward in the literature has been that the 
aforementioned cosmic strings could be the source of this 
CGWB [2l[. Besides CMB anisotropies, cosmic strings 
would also be (active) sources of gravity waves. There 
have been several works trying to obtain the GW spec- 
tra coming from cosmic strings, although there is still 
some debate in the community [2214281 ]. A full field the- 
ory simulation would be very helpful. The loop distribu- 
tions of strings is not completely understood either [29i — 
l39l ]. and that parameter could be key in estimating the 
gravity wave production from strings. Besides, the decay 
products of strings are also important, since strings could 
be decaying directly into gravity modes, or could decay 
into other particles. For example, pulsar timing bounds 
on cosmic string tension can change significantly due to 
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those factors mentioned above, i.e., loop production and 
decay mechanisms. 

All in all, the excess energy stored in relativistic com- 
ponents could come from gravity waves originated in cos- 
mic strings. In [2lj |. the authors discussed this possibility, 
and gave an upper bound on the cosmic string tension 
G^i < 2 x 1(T 7 at 95% confidence level (C.L.). How- 
ever, that zero order approach can be improved by not- 
ing the following: if the excess in N e g comes from cosmic 
strings, cosmic strings were there in the early universe. 
Then, cosmic strings would have also contributed to the 
temperature (and polarization) CMB anisotropics, and 
therefore the parameter estimation should be done in- 
cluding the contribution from string from the beginning. 

It is known that cosmic strings have degeneracies and 
correlations with other parameters in the ACDM model. 
For example, it was shown in that a n s = 1 was pos- 
sible with the (then) current data if strings were included 
in the analysis, due to a suit of degeneracies. Thus, it 
is interesting to investigate whether such correlations ex- 
ists between the string contribution and N e g, and what 
implications this may have both in the string tension and 
the extra neutrino species. 

This is the task we undertake in this paper: we per- 
form a Markov chain Monte Carlo type analysis to es- 
timate cosmological parameters, including cosmic string 
contributions, letting N e g be a free parameter. We use 
CMB data both for relative small I (WMAP7 [H) and 
for larger I (SPT [r|), as well as non-CMB data (BAO 
and HST HH). In section [TT] we describe the basic 
ingredients for our analysis, then we present our results 
in section UTTl and then conclude in IIVI 



II. METHODOLOGY 

In this section we discuss what information the number 
of effective neutrinos parameter A ff carries, then intro- 
duce the cosmic string template, and lastly, we show how 
the data analysis is done. 



A. Effective neutrino number 

The amplitude and the position of the acoustic peaks of 
the CMB power spectrum are determined by the energy 
content of the Universe before recombination, i.e the mat- 
ter and radiation content. So in this context it is crucial 
to have a refined knowledge of each of the components 
that could contribute to the radiation energy density at 
that epoch. The main component are photons, f2 7 , but 
other relativistic species like neutrinos can contribute to 
the radiation in the form: 
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where any extra radiation component is encapsulated 
in N Q ff, usually known as number of effective neutri- 
nos. One of its main effects is the change of the 



matter-radiation equality redshift, 



-a i 



which increases 



(decreases) as N c g gets smaller (bigger) . The same effect 
occurs if the energy density of matter, Q m h 2 , increases 
(decreases) . The parameter N e g can be expressed as [43| : 



N eS = 3.04 + 7.44 



n m h 2 3139 
0.13081 



- 1 



(2) 



Both those quantities, z eq and fl m h 2 , can be obtained 
directly from the matter power spectrum, allowing us to 
constraint N e g. 

The high energy physics standard model predicts three 
neutrino flavours, giving N c g = 3.046 which also ac- 
counts for some corrections from incomplete decoupling 
[44| . However, analysis on the N g from recent cosmo- 
logical data gives values greater than 3.046 indicating 
the presence of an extra relativistic energy component at 



early universe: a value of N e g = 4.34± 



+0.86 
-0.88 



at 68% C.L. 



for the combination of CMB data from WMAP— 7year 
with BAO presented in [43[ and the value of Hq given 
by the Hubble Space Telescope (HST) 0. Also the 
inclusion of CMB data at small scales, which provides 
tighter constraints, still shows a preference for the exis- 
tence of an extra component of radiation different from 
neutrinos. In combination with WMAP7 data, BAO 
and HST, an N cS = 4.56 ± 1.5 at 68% C.L is measured 
from the Atacama Cosmology Telescope (ACT) [l|| and 
iV eff = 3.86 ±0.84 with South Pole Telescope (SPT) data 

a 

This excess of relativistic energy could be an indication 
of an extra neutrino, but it could also indicate that there 
were other types of relativistic components in the early 
universe: it has been pointed in [19( this excess can be 
explained with the contribution of the primordial CGWB 
to the total energy density of radiation. 

As detailed in Sec. 35.7 of (2(j, it is possible to ob- 
tain an effective energy momentum tensor for the gravi- 
tational waves as an average of the stress energy carried 
by several wavelengths 



T, 



gw _ (hjk^hjk^) 
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where represent the covariant derivative of the ten- 
sor perturbation of the metric, hjk- For perturbations 
inside the horizon, gravitational waves can be considered 
to be propagating in a flat, Minkowski background. In 
this case, the equation of motion for the tensor perturba- 
tion of the metric becomes, d° 'd^h^ = 0, whose solution 
is a plane-wave = TZ{A^e lkaX }, with a wave vector 
k. This allows us to obtain an expression of the effective 
energy-momentum tensor: 



rpGW 



32tt 



(4) 
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which is the same as the energy-momentum tensor of a 
beam of non-interacting massless particles. Thus, the 
effects on the CMB and the matter power spectra of 
CGWB are equal to those produced by massless neu- 
trinos. This allows us to translate any constraint on the 
N e s into an upper bound on the gravitational wave en- 
ergy density, fl gw h 2 for those with a frequency larger 
than ~ 10 _15 Hz which corresponds to size of the sound 
horizon at decoupling. We can then obtain an upper limit 
of the GW energy density, 17 gw /t 2 , integrated over all the 
possible frequencies, see [45( for further details: 

f d (In f)h 2 Q gw (f) = 5.6 x l()- 6 (N cS - 3.046). (5) 
Jo 

The study of this extra radiation component is of 
special interest, since it carries exclusive information 
about the state of the primordial Universe, providing and 
unique window to explore the evolution of the Universe at 
those times, that will never be accessible by other ways. 
But, even if it is gravity waves, where do they come from? 



B. Cosmic string power spectra 

Cosmic strings are one type of cosmic topological de- 
fect that could be formed in the early universe; and if 
formed, they would still be around us. In high-energy 
physics models of inflation, cosmic strings arc ubiquitous, 
since they appear generically in models of SUSY grand 
unification and in brane inflationary models. 

The imprint of cosmic strings in the CMB is a diffi- 
cult problem to solve, because cosmic strings are highly 
non-linear objects that are active sources of the pertur- 
bation, so they keep sourcing the anisotropics; unlike in- 
flationary contributions, which act only around decou- 
pling and propagate freely since then. The thickness of 
cosmic strings is microscopical, though their length can 
be cosmological, and the time in which they evolve is 
also cosmological. It is thus very hard to get an analyti- 
cal handle on them (though there are works on effective 
theories of defects, e.g., |46|-[49|]) and numerically is also 
very difficult. There are different approaches to numeri- 
cally simulate cosmic strings: one can use the aforemen- 
tioned effective theory (known as the "unconnected seg- 
ment model", e.g., H [jj, Wtfj$ ) , one can forget about 
the microphysical details of the core of the strings and 
treat it as an infinitely thin object (the so-called " Nambu- 
Goto model", e.g., [ll M, M, M, HI ) , or one could try 
to perform a full field theoretical simulations to obtain 
the CMB power spectra ( which is usually done using 
an Abelian Higgs type model, e.g. [551 - 1601 ] ) . All those 
approaches have advantages and drawbacks, but in some 
sense they all are complementary. In this work, we have 
chosen to use the Abelian-Higgs type approach, and the 
string power spectra used is the one presented in [HD, [5!| . 

In order to perform the parameter estimation including 
cosmic strings, one would need to calculate the spectra 



of the strings at every point in the Markov chain. That is 
computationally very costly, so instead the cosmic string 
spectrum is computed one for a given parameter set (see 
[59( ) , and then the normalization of the spectra is left as 
a free parameter. This is justified in previous works, but 
it is easily seen as the contribution from the strings is 
subdominant. 

The normalization is given by the parameter G/i, where 
G is Newton's constant, and fi is the string tension. 
There is an analogous parameter that has been used in 
the literature previously: fio- This is the fraction of the 
total power coming from strings at multipole I = 10. Al- 
though the exact correspondence between Gfi and fio is 
model dependent, it is worth noting that /io oc (Gfi) 2 . 

The contribution from cosmic strings to the CGWB 
is even less understood than their CMB imprints. The 
production of GW from strin gs h as been studied for a 
number of years now [HI, 1231 l26l. l6l| . where the cusps 
and kinks on the strings where considered, mostly in a 
Nambu-Goto approximation. For those predictions, the 
size and distributions of loops are also of crucial impor- 
tance. One of the main focus of some of those works was 
to study how the GW contributions from strings would 
change when considering cosmic supcrstrings, and more 
specifically the probability of intercommutation of the 
strings p. For the present work, we only consider soli- 
tonic cosmic strings, and will consider p = 1, that is, the 
string segments always reconnect when they meet, which 
is an excellent approximation 1 . 

Even though the form of the GW spectra coming from 
strings is still under study (see (65;] for a different ap- 
proach), as is also the size and distribution of loops 
[36l 1371 log , for the present work we will use the ap- 
proach put forward in [27|, H8[ since it is one of the few 
that gives a direct translation from Gfi to gravity waves 
(and hence, to N c g) . Those authors show that under 
adiabatic initial conditions, incoherent superposition of 
cusp bursts from a network of cosmic strings (or supcr- 
strings) can produce a spectrum of gravitational waves. 
Assuming a certain model, it is possible to convert any 
bounds on Q &w h 2 into another over the string tension, 
Gfi [13, Us}; or v i ce versa. In [ll| the authors follow 
the procedure given in pjj which provides an analyti- 
cal approximation of the CGWB spectrum produced by 
strings. This approximate model allows us to relate the 
constraints on the Gfi to the gravitational wave energy 
density ( the string reconnection probability p will be 
fixed to 1 for our study). In [Hj] two different expression 
for the CGWB energy density are given depending on the 
string loops size. For those smaller than the horizon we 
have: 



This is the usual assumption, which is an approximation, albeit 
an extremely good one. There are however situations in which 
the probability of intercommutation of solitonic strings is not 
one llUH 
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Ogw(/) w 5 x 10" 2 G/t 



(6) 



whereas the spectra for the CGWB energy density given 
by loops of the size of the horizon is 



^gw(/) « 3.2 x VOT^y/Gn. (7) 

Therefore, once we have a value for G/i we can link it to 
a il gw h 2 which in turn is linked to an effective number of 
rclativistic species Nq^ whose effect would be analogous 
to the one produced by the gravity waves produced for 
those strings, all according to (28[. For subhorizon string 
loops the relations becomes 



G/xln 



pi 



5.6 x 10~ 6 N§^ 



5 x 10- 2 /i 2 

and for horizon-sized loops we have 
(G/x) 2 M pl 



G/i In 



3.6 x 10" 18 Hz 



5.6 x l(T 6 jVg^ 
3.2 x 10- 4 /i 2 ' 



(8) 



(9) 



In [2l|, Sendra and Smith calculated a best fit val- 
ues of N e g obtained in the parameter fitting to obtain 
a corresponding bound in G/i using the formulas above. 
In this work we follow the reverse path: starting from 
the values obtained of G/i, we use the formulas above to 
estimate the amount of gravity waves the strings would 
produce, and the effective number of relativistic species 
they would correspond to. 

Note that in this work we use several different effective 
number of relativistic species. We have N e s, which is the 
usual number used in the literature: this accounts for the 
effective number of any relativistic species at Big Bang 
Nucleosynthesis, except for photons and electrons. As we 
just mentioned, the cosmic strings present in our models 
will contribute contribute to the total relativistic species. 
The effective number of relativistic species analogous to 
the gravity waves produced from strings is what we call 
Nqw- This contribution is one of the contributions in- 
cluded in N c g. It is necessary to stress that there is a 
very high level of uncertainty in this number; not only 
because it will differ between long or short loops, but 
mainly because the translation from G/i to is far 

from understood. Lastly, since high energy physics tells 
us that there are three families of neutrinos, the abso- 
lute minimum number of extra relativistic species in the 
early universe must be 3.046. Thus, our models have two 
different sources for the effective number of relativistic 
species: the three neutrinos and the strings. Therefore, 
a new parameter is defined: N e g * = N c s — 3.046 — 
which accounts for any other extra ingredient necessary 
to fit the data best. A negative N e g* would tell us that 
either less than three neutrinos are preferred, or (more 
likely) that contribution from cosmic strings is too high. 



Parameter 


Symbol 


Prior 


Baryon density 


Q b h 2 


[0.005, 0.1] 


Cold Dark Matter density 


tl c h 2 


[0.01, 0.99] 


Angular size of sound horizon 


e 


[0.5, 10] 


Optical depth to reionization 


K 


[0.01, 0.8] 


Scalar spectral index 


n s 


[0.5, 1.5] 


Amplitude of scalar spectrum 


ln(10 1( M s ) 


[2.7, 4] 


Effective number of Neutrinos 


AW 


[1.047, 10.0] 


Cosmic Strings normalization 


(G M ) 2 


[0, 4 x 10" 12 ] 


Poisson point source power 


-^3000 


[0, 100] 


Clustered point source power 


n CL 

-^3000 


[0, 100] 


SZ power 


n sz 

-^3000 


[0, 100] 



Table I: Priors for the parameters used in our analysis, in 
the most general case with all parameters left free. For the 
different models we considered, some of this parameters will 
be kept fixed, as explained in the text. 



C. Data analysis 

The model considered is that of a ACDM with the 
usual six parameters: Amplitude of scalar spectrum A s ; 
Baryon density Cold dark matter density Q, c ; Optical 
depth to reionization k; Angular size of sound horizon 9; 
and the Scalar spectral index n s . We will refer to this 
model as the Power Law (PL) model. In our analysis, one 
or two more parameters will be added to the PL model: 
the effective number of neutrino species N e f{ (or equiv- 
alently the extra radiation component) and/or cosmic 
strings (parametrized as (G/i) 2 ). As we will show later, 
for certain cases we find that the data prefers n s = 1, in 
a similar way to what the authors in 0, 0, 0, 167} found. 
Therefore, we will also consider a model with n s fixed 
to one, and thus having only five free parameters. We 
will refer to this model as the Harrison- Zcl'dovich (HZ) 
model. Once again, N e s and G/i will be added to the HZ 
model. The linear priors used for the parameters can be 
found in t able [I] 

Besides the parameters mentioned above, on the fol- 
lowing tables we also show the values of some derived 
parameters: the Hubble parameter h; the contributions 
of strings to the total power spectrum at I = 10 /io; and 
two extra relativistic numbers, Nq^, which is the con- 
tribution to the effective number of extra species solely 
coming from strings, and N e g* , which is the contribution 
needed once the three neutrinos and the cosmic strings 
contribution has been subtracted. 

Our models will be compared with both CMB and 
non-CMB experimental data. Specifically, the CMB data 
we will use will be that of WMAP7 gfl| and SPT [H . 
When data from the SPT experiment is used, foreground 
contaminants have to be taken into account, therefore 
three extra parameters are used (normalized, as usual, at 



I = 3000): Poisson point source power D|oocb Clustered 
point source power D^ q \ and the Sunyaev-Zel'dovich 
power -Dfo^o- The non-CMB experiment data that we 
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will consider are the HST measurement of the Hubble 



constant (Hp) by Riess et al [42[ and the BAO by Perci- 
val et al I41L 



III. RESULTS 

The results of fitting WMAP7 data to the PL model 
with N e g and/or Gfi is shown in table ITU When we only 
allow N e g to be a free parameter and set G\x to zero, it 
is clear that, as it was pointed in previous works [l7lfl8l|. 
the preferred value for N c g is larger than 3, pointing 
out that there is a preference for some extra radiation 
component in the early universe. When including cos- 
mic strings also into the analysis, the preferred value is 
found to be still larger than 3, but smaller than when 
considering only N c g. Thus, that not only would cosmic 
strings be the source for the extra radiation component 
by means of creating gravity waves, but they also lower 
the need of the extra component, that is, G[i and N e g are 
anti correlated. For comparison, we give the numbers for 
the case where N e g is fixed to 3.046, and only the string 
contribution is allowed to vary. The likelihood for the 
models with 7 parameters (that is, the models that only 
have a string component, or only have a N e g component) 
are very similar, with a slight preference for the model 
with N e g. 

The value of N e g is quite poorly constraint by using 
only data from WMAP7 so other data are included. 
As mentioned before, we use HST and BAO data, and 
the results can be seen in tableHU The BAO and the HST 
data push most of the parameters in different ways; more 
specifically, BAO pushes the string contribution down 
and the N e g up; whereas HST pushes N e g down and the 
string contribution slightly up. When considering both 
data together, the tension between the two data comes 
into play, but the contribution from HST is more im- 
portant in that overall value of the parameters are tilted 
towards the preferences of HST. Comparing the result 
from fitting to WMAP7 only, to WMAP7+HST+BAO, 
we see that the mean value for the extra neutrino species 
is significantly lower. 

The tightest constraints on the number of effective neu- 
trinos so far is given by the inclusion of high I (SPT data 
P^D- The results can be found in table Hill both with 
and without non-CMB data. It is noticeable that when 
adding the high I data to the analysis, the effective num- 
ber of neutrinos decreases, similarly to what happened 
when non-CMB data are included. Another difference 
is that the string contribution is favored in this case: 
when comparing the model with only a string contribu- 
tion (PL+G/i), and the model with only N e g (PL+N e g), 
the likelihood is better for the model with strings (note 
that both models have the same number of parameters). 
Moreover, the model that fits both strings and neutrinos 
gives a 2a detection for strings. 

One striking difference from the case where only 
WMAP7 data was used is the fact that now the string 



contribution and the effective number of neutrinos are 
correlated, that is, when the contribution of strings is 
pushed up, so is the effective number of neutrinos. Fig- 
ure[T]shows the 2D likelihood plots for some combinations 
of the parameters, for two sets of data: on the one hand 
WMAP7+HST+BAO, and on the other WMAP7+SPT. 
The figure shows clearly that whereas fio and U^h 2 are 
correlated in the same way for both datasets, N e g is cor- 
related or anti correlated with the other two, depending 
on the datasets that are being fitted for. Thus, when 
fitting for the WMAP7+SPT dataset, an increase in the 
number of neutrinos brings an increase in the string con- 
tribution. Likewise, by allowing for string to be present 
in the analysis, we not only pick up a string contribu- 
tion, but the string contribution prompts the neutrino 
contribution to rise. 
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Figure 1: Two dimensional likelihood contours, where de- 
generacies between /k> and N e g and f2b/i 2 when fitting for 
WMAP7+HST+B AO , and when fitting for WMAP7+SPT. 
The shaded contours correspond to 1-cr and 2 — a contours. 
The figure shows that ft^h 2 and fio keep the same correlation 
for both datasets; but N e ff is correlated with the other two 
for one case, and anti correlated for the other. 



This trend can also be seen when including non-CMB 
data. In the case studied earlier, when we were fitting for 
WMAP7 data, the inclusion of BAO pushed the strings 
down, since it preferred to have a high N e g contribution. 
In the present case, when we are fitting for WMAP7 and 
SPT data, this preference for higher N e g means that the 
string contribution also goes up. In any case, as before, 
when both HST and BAO data are considered, the pa- 
rameters tend to favor the HST preferences. 

The inclusion of the non-CMB data, and most impor- 
tantly of the HST prior, constraints the parameters dras- 
tically better, and most of the degeneracies are broken. 
Figure [2] shows two dimensional likelihood contours for 
some of the parameters, when fitting for WMAP7+SPT 
with and without non-CMB data included. It is very 
apparent how the available parameter space has shrunk 
considerable, and also how the degeneracies are broken, 
specially between N e g and /iq. Note that the 2-a "dctcc- 
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Data 


WMAP7 


+HST 


+BAO 


+HST+BAO 


Model 


PL+G/i 


PL + iV cff 


PL+G/i+iV eff 


PL+G/i+AW 


lOOOb/i 2 


2.4±0.1 


2.21±0.06 


2.3±0.1 


2.32±0.09 


2.29 ± 0.07 


2.31±0.01 


Q. c h 2 


0.11±0.06 


0.18±0.03 


0.15±0.04 


0.118±0.02 


0.19 ±0.03 


0.12± 0.02 


6 


1.041±0.003 


1.029±0.005 


1.034±0.007 


1.039±0.006 


1.028± 0.004 


1.038± 0.006 


K 


0.10±0.02 


0.09±0.02 


0.09±0.02 


0.09±0.02 


0.09 ± 0.02 


0.09± 0.02 


n s 


0.99±0.02 


1.00±0.02 


1.00±0.02 


0.99±0.02 


1.02 ±0.02 


0.99± 0.02 


ln(10 1( M s ) 


3.18±0.05 


3.12±0.04 


3.10±0.07 


3.12±0.06 


3.13 ±0.05 


3.13± 0.07 


10 12 (G/x) 2 


0.18 (< 0.37) 


- 


0.15 (< 0.33) 


0.15 (< 0.33) 


0.10 (< 0.26) 


0.15 (< 0.34) 


N eS 


- 


7±2 


6±2 


3.6±0.9 


8±2 


4± 1 


h 


0.74±0.04 


0.84±0.08 


0.83±0.08 


0.74±0.02 


0.88 ± 0.07 


0.75± 0.022 


ho 


0.05 (< 0.107) 




0.04 (< 0.096) 


0.03 (< 0.091) 


0.03 (< 0.073) 


0.04 (< 0.096) 


N cs 


0.29/2.26 




0.26/2.16 


0.26/2.16 


0.22/1.94 


0.26/2.16 


N eS * 


-0.29/ - 2.26 


3.95 


2.69/0.79 


0.29/ - 1.61 


4.74/3.01 


0.69/ - 1.21 


Likelihood 


3752.780 


3752.619 


3752.475 









Table II: Marginalized likelihood constraints on model parameters, for three different models and with different datasets. 
The models differ in that N e fi and /or Gfi are parameters that we fit for, or are fixed parameters. The first three columns 
corresponds to the fitting to WMAP7 data, whereas the next three correspond to WAMP7+HST, WMAP7+BAO, and 
WMAP7+HST±BAO, respectively. The table shows the 6 usual parameters, plus the cosmic string contribution Gfi and 
the extra radiation component N e g. We also give the derived parameters h (the Hubble parameter) and /io (the contributions 
of strings to the total power spectrum at i — 10). Following the discussion in III Al and IHBI the values of and N e g* are 

given in pairs: the first number corresponds to subhorizon string loops and the second to horizon size loops. The values shown 
are the means and standard deviations, whereas the likelihood shown corresponds to the case that best fits the data. In the 
cases only upper limits can be placed, we give the mean value and next to it the 95% confidence limit in parenthesis. We do not 
show the likelihood of the three last columns, since the datasets used are different and there is no way of directly comparing 
the goodness of fit among them by only the likelihood. 



Data 


WMAP7±SPT 


±HST 


+BAO 


±HST±BAO 


Model 


PL±G/i 


PL ± iV cff 


all 


PL±G/x±AW 


lOOQb/i 2 


2.21 ±0.05 


2.27±0.05 


2.37±0.08 


2.25 ±0.05 


2.40±0.07 


2.25±0.05 


fi c h 2 


0.108 ±0.005 


0.13±0.01 


0.15±0.02 


0.13 ±0.01 


0.14± 0.01 


0.13±0.01 


6 


1.040 ± 0.002 


1.040±0.002 


1.036±0.002 


1.038 ± 0.002 


1.037±0.002 


1.038±0.002 


K 


0.08 ±0.01 


0.09±0.02 


0.09±0.02 


0.09 ± 0.02 


0.09±0.02 


0.08±0.01 


n B 


0.96 ±0.01 


0.98±0.02 


1.01±0.02 


0.97 ±0.01 


1.01±0.02 


0.98±0.01 


ln(10 10 A s ) 


3.17 ±0.05 


3.14±0.07 


3.06±0.07 


3.16 ±0.04 


3.01±0.06 


3.16±0.04 


10 12 (G M ) 2 


0.11 (< 0.22) 




0.24±0.09 


0.14 (< 0.27) 


0.26±0.06 


0.14 (< 0.25) 


N eS 




3.9±0.6 


6±1 


3.9 ±0.5 


5.7±0.9 


4.0±0.5 


n sz 

^3000 


6±3 


6±3 


5±3 


6±2 


5±3 


6±3 


n PS 
^3000 


20 ±3 


20±3 


19±3 


20 ±3 


21±3 


20± 3 


n CL 


5 ± 2 


5±2 


5±2 


4±2 


6±2 


5±2 


h 


0.72 ± 0.03 


0.75±0.04 


0.87±0.07 


0.75 ± 0.02 


0.89±0.06 


0.75±0.02 


ho 


0.03(< 0.057) 




0.07±0.03 


0.04 (< 0.072) 


0.08±0.02 


0.04 (< 0.067) 


N cs 


0.23/1.99 




0.34/2.44 


0.26/2.12 


0.35/2.49 


0.26/2.12 


N eB * 


-0.23/ - 1.99 


0.85 


2.61/0.51 


0.59/ - 1.27 


2.30/0.17 


0.69/ - 1.17 


Likelihood 


3763.980 


3764.717 


3760.712 









Table III: Analogous table to Table ILT1 but in this case the CMB data sets used are WMAP7 and SPT data. Because of the 
presence of the SPT data, we have to incorporate 3 more parameters, responsible for taking care of foreground effects. 



tion" for strings is lost when non-CMB data are included. 

The shrinking of the parameter space can be under- 
stood if one considers what is happening to the Hubble 
parameter. The inclusion of HST data basically narrows 



the prior space for the Hubble parameter Hq ~ 74.2 ±3.6 
km s _1 Mpc -1 , and the values of h that the parameter 
fitting prefers for the cases where the HST data are not 
included are much higher. Figure [3] shows the 2-D like- 
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Figure 2: Two dimensional likelihood contours for some of 
the parameters in the model. The shaded regions correspond 
to 1 — a and 2 — a confidence levels, where only CMB data 
are used for the analysis (i.e., WMAP7+SPT). The smaller 
contours depicted by solid and dashed lines correspond also 
to 1 — o and 2 — a confidence levels, but when also non-CMB 
(i.e., BAO+HST) data are included.. 



lihood contours for h versus N e g and /io- The shaded 
lines correspond to only CMB data, whereas the solid and 
dashed lines take into account also the HST+BAO data. 
It is clear that the HST just catches the lower end of the 
contour, and it is this effect which shrinks so noticeably 
the parameter space. 




70 80 90 70 80 90 



H o H o 

Figure 3: Two dimensional likelihood contours for iVeff and 
/io with respect to Ho. The shaded region represents the la 
and 2cr result obtain when only CMB data are used. Note that 
the mean of -Ho is higher than the value that HST experiments 
give (Ho ~ 74). That is why the allowed parameter space 
shrinks considerable when HST data are also included (solid 
and dashed lines). 

Figure[5]also shows that high vales of n s are preferred. 
This is important for inflation models coming from hy- 
brid supersymmetric inflation which seems to prefer val- 
ues of n s larger that the ones obtained in usual analy- 
ses. In fact, n s = 1 lies right in the center of the con- 
fidence contours for the CMB only case, and inside the 
2 — a level for the other case. A similar phenomenon 
was studied in some previous papers [H, 0, 0, where 
the Harrison-Zcl'dovich model was then considered as a 
viable model to explain the data. We have performed 
the same exercise, and the result is shown in Table IIV1 
where a model with n s = 1 (HZ) plus strings and N e ff is 
used to fit the CMB data (WMAP7+SPT). The values 
of the parameters change slightly from the case where n s 
was a free parameter, all within 1— a, which is expected 
since n s = 1 was right in the middle of the confidence 



Data 


SPT 


Model 


HZ+G^i+iVcff 


lOOfib/i 2 


2.31±0.04 


n c h 2 

e 


0.14±0.01 
1.037±0.002 


K 


0.09±0.01 


n s 


1 


ln(10 10 A s ) 


3.11±0.04 


10 12 (G>) 2 


0.19±0.07 


iVefi 


5.3±0.6 


n sz 

-^3000 


6±2 


^3000 


20±3 


n CL 
±73000 


5±2 


h 


0.84±0.02 


fw 


0.05±0.02 


/V cs 
' V GW 


0.29/2.30 


N eB * 


1.96/-0.04 


Likelihood 


3760.837 



Table IV: Analogous table to tables [TT1 and Hill In this case, 
the model we are fitting for is a HZ model, i.e., the scalar 
index is fixed to one n s = 1. The data we are fitting for is 
CMB only (WMAP7+SPT). This model fits the data better 
than any of the other models considered so far, even taking 
into account the number of parameters. Including non-CMB 
data would render this model less successful, partly because 
of the high value of Ho that the model prefers. 



contours. Thus, it is not surprising that in this case 
the value of the likelihood is very similar; but bear in 
mind that the number of parameters is different: PL has 
one more parameter than HZ. If one wishes to compare 
models with the same number of parameters, one should 
compare HZ+G/.i+N g with PL+G/i or with PL+iVeff; 
and in both those cases, the likelihood of the former is 
much better. Thus, the 10 parameter model (including 
the 3 parameters from SPT) that fits the data best, is 
the model that has n s = 1 and includes N e g and strings. 
Note, however, that also in this case the value of h is 
rather higher than what HST wants, so the inclusion of 
non-CMB data will disfavour this model. 

In all the analysis above there is one point of view 
that has not been taken into account: the strings that 
the different models prefer do in fact produce gravity 
waves which will contribute to the N e g, but that extra 
contribution is not taken into account as yet. In all three 
tables we have included a line giving Nq^, which is the 
effective number of relativistic species analogous to the 
gravity waves coming from the strings. The two num- 
bers that we are quoting are the ones obtained following 
the procedure of [2ll . |28| for subhorizon sized loops and 
horizon sized loops, and the big difference in both num- 
bers is already apparent. Note that we have chosen the 
mean value of (G/j,) 2 to obtain this numbers. Moreover, 
as we have mentioned before, this is a very rough esti- 
mate; therefore wc only use it as a hint of what could be 
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happening and the rich information we could obtain with 
it once the uncertainties are more under control. 

In other words, the N e g that we are obtaining from the 
parameter fitting has three possible sources: one that 
comes from the three species of neutrinos, another one 
given by -^VqW' an< ^ a third one that is not account for 
in the ingredients we have input in our model. The pa- 
rameter N c g* gives a measure of the level of those un- 
accounted effects. Most cases show the need of yet a 
contribution from another relativistic particle or source. 
Most notably, the predictions for smaller sized loops need 
a higher level of N e g*. In some cases N c g* is quite low; 
if it were exactly zero it would mean that the model is 
able to produce all the N e g it needs to fit the data. On 
the other hand, a negative N e f[* is an indication that the 
model is not good: the amount of cosmic strings pre- 
dicted is too high since their gravity wave production is 
too high. 



IV. DISCUSSION 

In this paper we have studied the correlation between 
cosmic strings and extra relativistic species (or effective 
number of neutrino species) by fitting different models 
to cosmological data. We considered the possibility of 
cosmic strings being the source of a cosmological gravity 
wave background, which in turn act as the extra relativis- 
tic species which the data seem to favor over the usual 
3 neutrino species. Cosmic strings are predicted in sev- 
eral inflationary scenarios, therefore, they seem perfect 
candidates to seed GW and account for the extra radia- 
tion. The idea of the cosmic strings being the sources of 
CGWB is not new, but this work expands the previous 
works in that the cosmic string contribution is included 
from the beginning, already in the parameter fitting pro- 
cess, and thus we are able to study different correlations 
between the different ingredients in these models. The 
inclusion of strings in the parameter fitting changes the 
value of N c g that the data needs. 

We have shown that the correlation between cosmic 
strings (G/z) and extra neutrino species (AT e ff) depends 
on the data sets used. When relatively low I data is 
used (WMAP7) these two components are anticorrelated; 
thus, the inclusion of cosmic strings into the model not 
only did account for the extra radiation species needed, 
but it actually lowered the need for them. However, when 
relatively higher I are included, the anticorrelation be- 
comes correlation, and an extra cosmic string compo- 
nent asks for a higher contribution from extra relativis- 
tic species. This (at first sight) unexpected effect may 
be understood by noting that a change in N c g changes 
the CMB temperature anisotropics power spectrum ba- 
sically in two ways: the height of the first peak and the 
position of the higher peaks. When only WMAP7 data 
are taken into account, only the height of the first peak 
is of importance, and cosmic strings are used to help in 
this endeavor. However, when higher £-s are taken into 



account, the position of the higher peaks gains more rela- 
tive importance, and this is counteracted by a correlation 
with cosmic strings. This corrclation/anticorrelation ef- 
fect not only happens between G\x and N e g, the same 
effect can be observed also between fit, and N e g. 

When non-CMB data are taken into account, the avail- 
able parameter space shrinks considerably. This effect 
comes mostly from the HST value for Ho. Typically, 
when the HST value of Hq is not taken into account in 
the analysis, rather high values of h are favored, and this 
is why adding the more stringent prior given by HST the 
parameters get much more constrained. The non-CMB 
data help constraint N e g considerably better when only 
WMAP7 data are used, and they break most of the de- 
generacies when also SPT data are used. The shifts in 
G[i and N e g due to these data sets depend also on the 
CMB data used, and on the correlation/anticorrelation 
effect. In general, HST tends to lower N e g whereas BAO 
tends to increase it. Therefore, when the only data used 
are the WMAP7 data, BAO tends to suppress the string 
contribution, but when also SPT data are used, tends to 
favor a higher string contribution. It would be interest- 
ing to analyze whether new CMB data alone, for example 
Planck data, could break the degeneracies between N e g 
and G/i without the need of non-CMB external data [68| . 

The analysis also shows several remarkable results. 
CMB data alone prefer a high value for the N c g, higher 
than the standard 3 neutrino species. However, when 
HST data are included, the results are consistent with 
only 3 neutrinos. In most cases, the string contribution 
is discrete, and no strings are consistent with the data ex- 
cept in the case where CMB+SPT are used with no-CMB 
data, where a 2a detection (or rather, hint) is obtained . 

Note also that in many cases an n s = 1 perfectly con- 
sistent with the data, in fact, we studied what is the 
goodness of fit for a HZ+G/i + N c g model, fitting to all 
the CMB data, and found that the model fits the data 
extremely well. Actually, the HZ +G/i + N e g model fits 
the data better than a PL +N c g or a PL+G/j, model , 
which have the same number of parameters. The like- 
lihood for a HZ +G/i + N c g is actually very similar to 
a PL+G^ + N Q ff , even though the latter has one more 
parameter. 

In order to fit the data a certain number of cosmic 
strings arc preferred, and those cosmic strings would cre- 
ate some gravity waves that are responsible for part of 
the N c g. We have tried to factor out this contribution, 
as well as the contribution of the standard three neutri- 
nos, to estimate whether some other source to N e g was 
still needed, and encoded it in N e g*. A value N e g* = 
would mean that all the players in the model are enough 
to account for all the N e g needed to fit the data, i.e., the 
three neutrinos and the cosmic strings present do the job. 
We find that in some cases the number can be somewhat 
small, but in many other cases the need for yet another 
extra relativistic source is needed. More importantly, 
sometimes N c g* is negative, hinting for some incongru- 
ence: the amount of strings chosen from the fitting is too 
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high. This happens in several of the situations analyzed 
by us; it could be a good measure to disqualify models as 
not viable, but unfortunately our feeling is that the un- 
certainties in obtaining are too high to make any 
kind of claim. One measure for the level of uncertainty 
could be the big difference in numbers for smaller and 
bigger cosmic strings loops; but all other uncertainties 
in the way those estimates are obtained may be much 
bigger. 

We should mention the shortcomings and possible im- 
provements for this work. The string power spectra used 
in this work is the one coming from Abelian Higgs field 
simulations. There are other approaches to obtain the 
string spectra, which have roughly the same form, but 
the details might change the results (though we may an- 
ticipate that only slightly). However, a better under- 
standing of the power spectra coming from strings are 
high £ would be good. As already mentioned several 
times, the main source of uncertainty comes from the 
production of gravity waves from cosmic strings. There 
is intense work in this subject from different groups, us- 
ing different approaches, and there is no consensus in the 
community. Unlike in the case of CMB predictions from 
strings, where different approaches give relatively fairly 
similar results, the loop production and decay process of 
strings is still a very open issue. In this work we have 
used one of the few works which give a recipe to translate 
from cosmic string Gfi to gravity waves, which in turn we 



transform into N c g, but we believe there are many un- 
certainties and assumptions that need to be checked an 
improved. We, therefore, consider this work as a first step 
into the analysis of the role of cosmic strings as sources 
of extra radiation component in the universe, which de- 
mands further understanding of the underlaying cosmic 
string dynamics. Moreover, it would also be interesting 
to perform an analysis including extended cosmological 
parameters as in |69(. 
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